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ABSTRACT 


The  importance  of  designing  a  radar  receiver  to  operate  effectively 
in  the  presence  of  noise  which  is  correlated  with  the  signal  creates  the 
necessity  for  developing  a  theoretical  model  of  the  clutter  which  can  be 
used  for  the  evaluation  of  various  detection  techniques. 

Kelley  and  Lerner  in  reference  1  have  developed  one  such  model. 
Reed,  in  reference  2,  has  developed  another  model.  This  report  presents 
a  refinement  of  Reed's  original  model  for  the  clutter  noise  process. 

The  clutter  cloud  is  considered  to  be  a  collection  of  individual 
point  scatterers  moving  acout  and  reflecting  signals  independently  of 
one  another.  The  individual  motions  of  the  scatterers  are  considered  to 
be  independent  of  one  anotner,  but  the  cloud  is  allowed  to  have  an  over¬ 
all  drift  velocity. 

The  analysis  technique  is  the  following:  An  RF  pulse  is  trans¬ 
mitted  at  t  =  0.  The  wavefront  is  reflected  from  the  scatterers  at  the 
distances  r^.  If  we  neglect  the  multiple  scattering  effect,  then  the 

returned  echo  consists  of  many  replicas  of  the  transmitted  waveform 

2r, 

beginning  at  times  t,  =  — £  ,  where  c  is  the  velocity  of  propagation 

K  C 

in  the  medium.  A  probability  density  that  a  scatterer  will  be  at  range 
r  and  time  t  is  assumed.  Since  the  scatterers  are  independent,  the 
arrival  of  an  echo  at  time  t  does  not  affect  the  probability  of  an 
echo  arrival  at  other  times;  thus  the  echo  arrivals  constitute  a  non¬ 
stationary  Poisson  process  with  rate  v(t). 

The  calculations  proceed  from  this  point  to  determine  the  non¬ 
stationary  expected  value  and  covariance  function  for  the  clutter  noise 


iii 


process. 

determine 


The  process  is  then  approximated  as  stationary  in  order  to 
the  power  spectrum. 
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SUMMARY 


The  time-varying  correlation  function  for  a  radar  (sonar)  back- 
scattered  noise  process  is  determined.  An  approximation  is  made  in  order 
to  obtain  an  expression  for  the  power  spectrum  of  the  noise  process.  The 
correlation  function  and  the  power  spectrum  are  both  a  function  of  the 
probability  density  that  a  scatter  will  be  at  range  r  and  time  t  as 
well  as  the  signal  pulse  shape.  This  implies  that  in  order  to  minimize 
the  noise  power  the  signal  may  not  be  wideband  but  is  affected  by  the  * 
probability  distribution  of  the  scatterers. 
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I.  INTRODUCTION 


In  order  that  a  radar  (or  sonar)  receiver  may  be  properly  designed 
to  operate  in  an  "optimum"  fashion  in  the  presence  of  clutter  (reverbera¬ 
tion)  it  is  necessary  to  determine  the  statistical  properties  for  the 
back-scattered  noise  process  (which  will  be  correlated  with  the  signal 
reflection  to  be  detected). 

The  statistical  properties  for  the  noise  process  that  are  analyzed 
in  this  report  are  the  time  varying  correlation  function  and  the  power 
spectrum  (for  the  steady  state  case).  Since  it  would  be  difficult  to 
determine  the  power  spectrum  for  a  noise-like  returned  echo  train  for  a 
particular  scattering  cloud,  the  problem  is  idealized  at  the  outset.  A 
refinement  to  Reed's  first  order  approximation  (reference  2)  to  the  power 
spectrum  of  the  returned  process  is  analyzed.  This  process  is  applicable 
for  a  wide  class  of  possible  scattering  clouds. 

The  results  of  the  present  analysis  are  an  extension  of  the  results 
obtained  by  Reed.  Reed  assumed  that  the  set  of  distributed  targets  is 
large  in  extent  and  is  composed  of  a  homogeneous  random  collection  of 
individual  scatterers  which  are  possibly  moving  about  independently  of 
one  another  with  possibly  an  overall  drift.  The  effects  of  multiple 
scattering  are  neglected  in  order  to  obtain  a  tractable  answer.  In 
addition  the  scatterers  are  assumed  to  have  small  velocities  so  that  an 
echo  from  a  particular  scatterer  can  be  regarded  as  a  doppler-shifted 
replica  of  the  transmitted  waveform.  The  result  Reed  obtained  for  the 
stationary  correlation  function  is  given  in  the  conclusions  of  this 
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report,  along  with  the  new  results  of  this  report.  A  comparison  of  the 
two  results  is  also  made  in  the  conclusion. 

The  next  section  discusses  the  assumptions  pertinent  to  the  analysis 
and  the  refinements  of  the  first  order  approximation  obtained  by  Reed 
before  determining  the  time-varying  correlation  function  of  the  noise 
process. 


II.  REFINEMENT  OF  THE  MODEL  AND  ITS  ANALYSIS 

The  following  assumptions  are  similar  to  Kelly  and  Lerner  in 
reference  1  as  is  the  method  of  calculation.  The  pattern  of  the  calcu¬ 
lations  follows  that  of  Reed  in  reference  2. 

While  most  of  the  following  assumptions  and  comments  have  been 
discussed  in  references  1  and  2  they  will  be  repeated  here  for  complete¬ 
ness. 

The  refinement  of  the  model  consists  of 

1.  the  inclusion  of  the  three  dimensional  aspects  of  the  problem 
as  well  as  the  statistical  structure  of  the  reverberation 
(clutter)  cloud, 

2.  a  coordinate  system  centered  at  the  transmitter  of  the  sonar 
(radar) , 

3.  the  inclusion  of  the  transmitter  antenna  g;  in  function, 

G(§,  0),  normalized  to  unity  in  the  direction  of  maximum 
gain  (fl  is  the  depression  from  the  vertical,  0  is 
azimuth) . 

The  effects  of  scatterer  rotation  as  well  as  multiple  scattering  will 
be  ignored. 
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Following  Reed  (reference  2),  if  the  transmitted  signal  is  given  by 
f(t)  *  s(t)  e2nlfo^  where  f(t)  is  the  real  transmitted  signal  plus  i  times  its 
Hilbert  transform  (reference  3i  p.  320),  then  the  returned  echo  from  a  collec¬ 
tion  of  scatterers  at  (r,  9,  0)  as  a  function  of  time  over  an  observation 
interval  (0,  T)  is  given  by 


G(9k,  0k> 


tke(0,  T) 


ak  s(t-tk)  exp  i[2nfk(t-tk)  +  6k  +  e(V  ^kO 


where 


0  includes  such  factors  as  RF  gain,  antenna  matching,  etc. 
a^  is  the  amplitude  reflection  coefficient 
,  2v  2vk. 

fk  •  ^  1  +  +  -jj— )  f0  is  the  doppler  frequency  (introduced  by 

Reed  in  reference  2) 


vq  is  the  drift  radial  velocity  of  the  system  of  scatterers 

vk  is  the  differential  radial  velocity 

6k  is  the  phase  shift  due  to  the  effects  of  reflection 

e(fl  ,  0  )  is  a  phase  shift  due  to  the  transmitting  antenna  beam 

K  K 

phase  characteristic 

c  is  the  velocity  of  sound  (or  light)  in  the  medium. 

The  subscript  k  refers  to  the  kth  scatterer. 


The  following  statistical  assumptions  concerning  the  random  values 
0k.  fk.  0k.  ak»  £>k  are  made: 

1.  The  random  values  are  taken  from  a  joint  probability  distribution 
where  it  is  presumed  that  a  scatterer  is  present  at  range 
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Following  Reed  we  next  assume  the  scatterers  are  Poisson  distributed 
in  range  (time).  It  has  already  been  noted  that  they  are  independent.  The 
joint  probability  distribution  for  the  Poisson  process  for  exactly  n 


events  to  occur  in  the  intervals  (t,,  t,  +  dt,) 
this  nonstationary  case  can  be  expressed  as: 


.(t  »  t  +  dt  )  for 
n  n  n 


(4) 


where  it  has  been  assumed  that  the  signal  duration  A  is  short  compared 
with  T,  the  observation  interval,  and  we  can  neglect  the  end  effects 
near  0  and  T. 


With  these  assumptions  we  may  now  determine  the  expected  value  of 


p(a)  da  J*.  ..J*  p(h) 


x  J...J  p(0)  d 0  J. .  .J  p(6)  d6  J...J  p(f)  df 


"  G(0.,  0  ) 

Z  - 2 -  “k  s(t_V  exp  i[2nfk(t-tk)  +  6R  +  eify  0R)] 

k=l  t  k 


where 


and  a 


-*  -•  n 

p(x)  dx  *  It  p(xt)  dXj^  , 


has  been  absorbed  into 


P. 
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If  we  assume  that  each  of  the  random  variables  is  stationary  except 
for  range  (or  delay  t)  or  vary  slowly  in  time  compared  to  the  pulse  repe- 


tition 

frequency 

,  then  E  [  z(  t)]  can 

be  expressed  as  follows: 

E  Tz(t)]  =  P  E(a 

00  ^ 

1  »i ,(1>I  I-J 7T 

n=0  0 

v(tk5  "4 

v(t)  dt 

dtk 

k 

x  j.:. 

•  oo 

!  pC  9)  de 

_  00  _  ^  _  00  _ 

p(0)  d0  J _ J 

•00  *00 

■*  w 

p(f)  df 

n 

V  V 

o(ek,  Kk) 

s(t-tk)  eir2n(t-tk) 

fk  ♦  E(9k* 

0J1 

x  L 

k=l 

K  K 

K 

(6) 


where 


(7) 


s  characteristic  function  of  6  with  u=l 


and  the  a,  's  are  identically  distributed. 

If  we  let  L(  9^ v  0^)  ■  G(sk*  0^)  e^e^k*  and  assume  that  the 

L's  are  identically  distributed  then 

r  pT  -  n  v(t  )  "JQ  v(t)  dt 

E[z(t)]  «  P  E(a)  E(L)  ^(1)  £  ^  — -j~*  e  dtR 

n*0  0  k=1 

x  J.T.J  p(f)  d7  l  e2niVt-tk)  (8) 

k-1  1  k 


The  approach  used  in  obtaining  equation  (8)  differs  from  that  used 
by  Kelly  and  Lerner  in  that  they  make  an  assumption  which  eliminates  the 
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doppler  shift  effect.  Consequently,  they  obtain  results  which  differ 
considerably  from  the  results  which  are  derived  in  the  following  pages. 


i-l  ,T 


dt'^  -J  v(t) 

E  [z(t)]  =  P  E(a)  E(L)  4^(1)  ^  - — - 


dt 


n  e 


n=0 


rT-it^)  { r  p(f) 

J0  t'  L  J 


e2nif(t-f)  df|  dt, 


dt'  "f1  -J  v(t) 

=  P  E(a)  E(L)  *6(1)  £  — - - 


dt 


n  e 


n=0 


x  f  4  •( t— t * )  dt' 

J0  (t')2 


=  P  E(a)  E(L)  <Ml)  f  s(t-t')  *(t-t' )  dt' 

6  J0(t')2 


where 


2nif(t) 


*(t)  -  J  p(f)  e^'"'”'  df 

»  characteristic  function  of  p(f). 


(9) 


(10) 


If  v(t)  is  replaced  by  n  y  t2  u  y  , 


then 


J  .3 


E[z(t)]  «  P  E(a)  E(L)  *6(1)  J  uCSF’lD  s(t-t')  *(t-f)  df  (11) 


where  u(r,t)  is  the  probability  density  of  finding  a  scatterer  at 

r3 

(r,  t).  Let  f(t-t')  -  s(t-t')  t(t-t'),  a  -  P  E(a)  E(L)  ♦  (1)  n  y 
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then 


provided  f(t-t')  -  f(t')  ~  0  (13) 

or  the  second  integral  is  approximately  zero. 

Equation  (13)  establishes  one  functional  relationship  for  f(t) 


which  may  be  useful  later  in  determining  an  optimum  signal  pulse  shape, 
s(t). 

If  the  interval  (0,  T)  is  large,  then  E[z(t)]  can  be 
approximated  by 

c* 

E[z(t) ]  ■  a  j*  ^  s(t)  #(t)  dt  (14) 

-•§ 

This  is  a  generalization  of  Cambell's  theorem  (see  references  1  and  2). 

In  a  similar  manner  we  can  show  that  the  covariance  function  of 

z(t)  is 

cov  £  z(t-t)  z(t)J  ■  E  £  z(t-T)  z(t)^  -  E  £  z(t-T)^  E  £  z(t)J 

T 

-  r  p !2  E(a2)  E(G2)  4(r)  f  s(t-T-t' )  s(t-t')  dt'  .  (15) 

J0  (t'T 
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9 
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However*  equation  (1?)  gives 


•  (£ 1  ts) 

n  ,  =  *  r  »Ll±i2  i  s(t)  i2dt . 

o  2  2  J 

•OB 

This  shows  that  the  noise  power  may  or  may  not  increase  linearly  wi  th 

T  .  since  the  probability  density  u(r,  t)  now  influences  the  noise 

s’ 

power.  Thus,  the  signal  pulse  shape  is  strongly  determined  by  the  pro¬ 
bability  density  u(r,  t)  if  the  noise  power  is  to  be  minimized;  i.e., 
for  one  type  of  probability  density,  w(r,  t),  the  signal  spectrum  might 
be  broad,  while  for  another  u(r,  t)  the  signal  spectrum  might  be  narrow 
It  is  shown  in  the  appendix,  in  a  similar  manner  as  that  used 
in  reference  1,  that  z(t)  approaches  a  Gaussian  process. 
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APPENDIX 


We  show  here  that  the  probability  distribution  of  the  returned 
echo  noise  waveform  approaches  a  Gaussian  probability  distribution  as  the 
mean  number  of  scatterers  per  unit  distance  becomes  large. 

*e  start  by  considering  Reed's  original  model  of  the  returned  echo 
noise  process,  which  is  a  simplification  of  equation  (1).  If  the  time  inter¬ 
val  of  observation  (0,  T)  is  large  compared  with  the  transmitted  pulse  dura¬ 
tion,  A,  a  returned  echo  waveform  in  (0,  T)  for  a  transmitted  pulse  is 


z(t)  =  y  sCt-t^)  exp  2nifk(t-tk) 
tk  e  (0,  T) 


where  the  various  symbols  have  been  defined  in  the  main  body  of  the  report. 
Following  Kelley  and  Lerner  [reference  11  we  can  write  immediately 


P(z)  =  \  f*  •  •  j  pCa)  da  [*•  •  *  [*  p(t)  d£  [*•  •*  p(f)  df  6(z-z  ) 

'  .  \S  J  .  •  Ii 


n=0  -« 


where 

n 

(1)  o(x)  d(?)  =  |l  ptx.)  dx.  =  joint  probability  distribution  of  the 

i  =  l  1  1 

x^'s,  where  x^,...,  xr  are  assumed 
independent. 

(2)  We  disregard  effects  within  A  of  the  end  points  of  (0,  T); 

n 

zn  =  y  a^  s(t-tk)  exp  2nifk(t-tk)  represents  the  waveform 

k-1 

resulting  from  events  at  the  times  t^  (0  <  t^  <  T); 


12 


z(t)  =  )  a^  s(t-tk)  exp  2nifk(t-tk). 

k=-« 

The  characteristic  function  of  the  distribution  is 

o» 

4>(u)  =  \  p(z)  U  Z  dz 

-® 

=  7  V  nJ  p(0  d(an)  M  d*n 

n=0  —oo  -on 

T  n 

.  f...j  dtn  exp  |^i  u  y  s(t-t^)  e2nifk^t-tk^~| 

0  k=l 

00  ^  00  T  00 

-  y  - — “J —  {  [  p(a)  da  f  dt’  [*  p(f)  exp£i  u  a  s(t-t')  e2nlf < t~t'  )"j  df], 

n=o  -oo  0  0 

T  00  a> 

=  expj-vT  +  w  J  dt'f  p(a)  da  J  p(  f )  df  expfiua  s(t-t')  e£,llf  ( t-t ' J  ~|  df! 

0  -00  —CO 

where  v  is  not  time  varying. 

If  we  expand  the  exponential  involving  s(t)  in  a  series  and  introduce  the 
characteristic  function  for  pC f )  as 

00 

0(t)  =  f  p(f)  e2nlft  dt 
—  00 

then 

00 

I(s,0)  =  f  p(f)  df  expfiua  s(t-t') 

—  00 

2 

=  1  +  iua  s(t-t')  0( t-t '  )  +  -  1y,a...r5^-.t '  ) 2  gfoct-t')] 
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Next  we  have 


f  p(a)  I(s,0)  da  =  1  +  y  - "~J~"  ar  sr(t-t')  0Fr(t-t')"j 


-ao  r=l 

r  r 

where  a  denotes  the  expected  value  of  a  . 

Using  the  above  we  obtain 

00  r  __  ^ 

(J)(u)  =  expjv  V  ar  J  sr(t-t')  0fr(t-t')l  dt'”)-  • 

r=l  0 

The  probability  density  of  z  is  found  by  transforming  its  charac¬ 
teristic  function  as  follows: 


Hz)  =  ^  j  <l>(u)  e  iuz  du 


00  oo 

=  ~  J  expj-iuz  +  v  Y  ar  J  sr(t-t' )  0f r ( t- 1 *  ) ~[  dt'1 


du 


r=l 


=  —  J  du  exp/-iu(z-z)|  exp-fv  ar  J  sr(  t-t '  )  0fr(  t-t' ) "]  dt'  | 


r=2 


=  ^  f  exp-j-iu(z-z)]-  expjv  Y  ~T“  ar  ar(t)T 


du 


r=2 


where 


3r(t)  ■  J1  sr(t-f)  r( t— t * )  1  dt'  . 


The  above  gives  the  time  varying  probability  distribution  of  z. 
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We  assume  next  that  the  noise  process  is  stationary  in  order  to  simplify  some 
of  the  following  calculations  and  be  able  to  express  the  probability  density 
of  z  in  terms  of  its  stationary  variance  and  mean.  Thus,  the  probability 
density  of  z  is  given  approximately  by 

P(z)  =  J  exp{-iu(  z-zl  exp{-~  expjv  ^  ^  du 

-od  r=3 

where  g  is  no  longer  a  function  of  time,  .ve  next  make  the  following  sub¬ 
stitutions 


z  -  z  =  \Z"2  x;  /~2  u  =  T 
z  z 


and  in  addition  we  recall  that  for  v  a  constant  then  V  2  =  v“  a.  Then 

z 


/"l"  =  v11  a. 


P(z)  dz  becomes 


P(x)  dx 


dx  T  f  fV  (i)r  Co  J  r  1 

*  IT  J  expilTl*  -  if  expl  L  TT  a  ar  \ 

-«  r=3  2 


dT 


If  we  expand  the  second  exponential  of  the  above  integrand  in  a  series  we 
obtain 

P(x)  =.  3  J  40  expf-iT*  -  {i  *  -  f,  ^  -JT1 


.£31 -Li  ,  CS 

v  +  5i 


ini  T 


LZ  Lis 

2 


•  } 
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Then  recalling  that  Tsee  reference  4,  p.  2281 


(n) (  v 
©  ( x; 


itx 


dx  =  (-it)“  e 


n=0,l,2, . 


(n) 

5> 


(x) 


dn  1 
dx11  ^ 


0(x)  =  0(o)(x)  = 


we  finally  obtain  the  Edgeworth  series 


P(x)  =  0(x) 


1 

K 


0(3)(x) 


TJ 

1  a  "4 

+  47  ~T~ 

a 


1 

v 


0(4)(X) 


Thus  P(x)  approaches  a  Gaussian  distribution  as  v  becomes  large,  where 
v  is  the  mean  number  of  events  per  unit  time.  If  p  is  the  mean  number  of 
scatterers  per  unit  distance  in  range  which  are  illuminated  by  the  transmitted 
pulse  and  intercepted  by  the  receiver,  then 


v 


Thus  the  returned  echo  noise  process  approaches  a  Gaussian  process  as  o, 
the  mean  number  of  scatterers  per  unit  distance,  becomes  large. 

It  is  a  straightforward  matter  to  use  the  above  formal  methods  to 
show  that  all  joint  distributions  are  Gaussian  in  the  limit. 
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For  the  case  when  z(t)  is  given  by  equation  (l)  and  v(t)  is  given 
by  equation  (3),  then  the  probability  distribution  will  again  be  Gaussian 
provided  the  comments  made  above  regarding  ar  apply  to  the  modified  ar 
given  by 

oT(  t)  =  ur  ^  sr(t-t' )  0^r(t-t' )  ~|  dt '  . 
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